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I.  onBjawa  of  fioibct  activitibs 

The  interection  between  e  vibrating  euboerged  structure  and  the  sur- 
rounding  fluid,  which  features  coupling  between  the  surface  pressure 
distribution  and  the  structural  displaceaent ,  is  an  inherent  feat\ire  for 
sound  radiation  and  target  strength  analyses.  A  variety  of  approaches  have 
been  iapleaented  in  the  past,  but  each  suffers  from  serious  liaitations. 
Fomal  BatheBatical  analysis  luing  separation  of  variables  or  integral 
transfom  techniques  is  suitable  only  for  the  siaplest  structural  aodels, 
while  full  finite  eleaent  descriptions  of  realistic  structures  and  the  sur¬ 
rounding  aediua  lead  to  excessively  large  coaputer  siaulations.  One 
approach  uses  approxiaate  iapedance-type  botmdary  condition  of  uncertain  ac¬ 
curacy  to  aodel  the  fluid  response.  Boundary  eleaent  foraulations 
rationally  represent  the  interaction  phenoaena  without  explicitly  solving 
field  equations  for  the  fluid,  at  the  expense  of  an  enomous  increase  in 
coaputational  effort  due  to  the  need  to  cover  the  surface  with  a  reasonably 
fine  aesh. 

Because  of  the  liaitation  of  classical  analytical  techniques,  it  was 
widely  believed  that  analytical -type  solutions  can  only  be  obtained  in 
idealized  systeas.  The  priaary  objective  of  the  project  was  to  develop  an 
analytical -type  approach  for  aodeling  fluid- structure  interaction  in  the 
frequency  doaain.  This  objective  was  achieved  by  developing  the  surface 
variational  acoustics  principle  (SVP) .  SVP  represents  the  spatial  dependence 
of  the  surface  response  as  a  series  expansion  in  a  set  of  assuaed  basis 
functions.  Since  the  solution  derived  froa  SVP  is  the  coefficients  of  this 
series,  the  priaary  difference  froa  the  results  of  a  classical  analysis  lies 
in  the  fact  that  SVP  deteraines  these  coefficients  by  nuaerical  techniques. 
In  essence,  iapleaentation  of  SVP  is  cMparable  to  using  the  aethod  of  as- 
sxsMd  Bodes  in  structural  dynaaica  When  the  aechanical  energy  functions  are 
very  coaplicated,  as  in  the  case  of  a  doubly  curved- shell,  deteraination  of 
the  systea  coefficients  for  the  latter  aethod  also  requires  nuaerical  tech¬ 
niques  to  evaluate  the  stiffness  and  inertia  coefficients.  A  specific 
advantage  of  SVP  foraulations  in  coaparison  to  conventional  finite  eleaent 
and  boundary  eleaent  .  alations  is  the  fact  that  convergence  of  the  rep¬ 
resentation  can  be  Jud;~d  by  exaaining  the  solution.  One  aerely  needs  to 
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Identify  that  the  higher  order  aodal  eaplitudes  ere  small  ••  there  is  no 
need  to  construct  e  more  refined  solution  to  use  as  a  reference  for  conver¬ 
gence.  In  ell  of  the  problems  addressed  by  SVP,  convergence  was  found  to  be 
quite  rapid.  Furthermore,  even  when  a  result  was  not  fully  convergent,  it 
was  found  nevertheless  to  be  remarkably  close  to  the  convergent  result. 

Because  SVP  represents  the  surface  pressure  in  a  functional  series 
form,  it  involves  a  substantially  reduced  set  of  unknowns  in  comparison  to 
bowdary  element  formulations,  which  ere  based  on  point-wise  discretiza¬ 
tions.  In  addition,  SVP  is  essentially  an  optimization  process  that  selects 
the  coefficients  of  the  modal  series  so  as  to  minimize  the  deviation  of  the 
derived  solution  relative  to  the  true  one. 

The  objective  of  the  first  phase  of  the  project  was  to  develop  the  SVP 
in  sittiations  where  the  surface  motion  is  specified.  Such  a  case  arises 
when  a  rigid  body  is  set  into  a  specified  motion.  Meeting  this  objective 
entailed  validating  SVP  results  against  known  analytical  solutions.  The 
next  research  phase  addressed  the  broad  objective  of  extending  SVP  to  treat 
acoustic  radiation  from  elastic  structures.  This  objective  was  achieved  by 
coupling  the  SVP  equations  to  the  method  of  assumed  modes  for  the  structural 
vibration.  The  latter  method  was  chosen  for  the  description  of  the  struc¬ 
tural  motion  because  it  too  is  variationally  based,  being  founded  on 
Hamilton's  principle.  In  the  first  two  research  phases,  the  analysis  was 
restricted  to  cases  of  axlsynmetrlc  excitation.  The  next  phase  met  the  ob¬ 
jective  of  generalizing  SVP  to  treat  situations  that  are  not  axisymmetric . 
Two  cases  were  considered:  nonsymmetric  excitation  of  bodies  of  revolution 
and  two-dimensional  (plane  strain)  models.  The  final  objective  was  to  ex¬ 
tend  SVP  to  treat  scattering  of  waves  at  arbitrary  incidence  on  an  elastic 
body. 


Two  considerations  influenced  the  selection  of  problems  for  each  re¬ 
search  phase.  Problem  that  had  been  solved  previously  provided  validation 
of  SVP.  Once  the  formulation  was  validated,  problems  were  selected  on  the 
basis  of  their  relevance  to  specific  structural  acoustic  phenomena.  A  par¬ 
ticular  problem  area  that  received  much  attention  was  the  influence  of 
eigenvalue  veering  phenomena  on  predictions  of  acoustical  performance.  Such 
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phenoMna  arise  when  alteration  of  a  systea  paraaeter  catises  the  natxiral 
frequencies  of  different  aodes  to  becoae  slallar.  The  in> vacuo  aode  shapes 
In  this  case  are  hl^ly  sensitive  to  the  value  of  the  systea  paraaeter, 
idilch  suggests  that  the  forced  response  of  the  subaerged  systea  should  dis¬ 
play  slallar  sensitivity. 


II.  BASIC  nmCIPLE 

The  derivation  of  the  variational  principle  is  fairly  stral^tforward. 
The  starting  point  Is  the  surface  Integral  equation  for  noraal  velocity. 

The  steps  Involved  In  deriving  this  eqtiatlon  begin  by  taking  the  gradient  at 
a  field  point  of  the  Kirchhoff-Helaholtz  Integral  theorea.  Before  the  field 
point  Is  brought  to  the  surface,  a  regularization  procedure  Is  used  to 
reaove  a  hypers Ingularlty.  This  regularization  Is  needed  becatise  one  would 
otherwise  need  to  Integrate  over  a  double  gradient  of  the  Green's  function. 
The  eqxiatlon  derived  In  this  Banner  is  called  the  second  surface  integral 
equation.  The  field  points  and  source  points  appear  In  synnetrlcal  fora  In 
this  Integral  equation,  so  It  Is  laalnently  suitable  for  constructing  a 
variational  principle.  The  SVP  theorea  Is  obtained  by  aultlplylng  the 
second  sxirface  Integral  eqxiatlon  by  a  virtual  Increaent  in  the  surface  pres¬ 
sure  at  a  surface  field  point,  and  then  Integrating  over  all  such  points  on 
the  sturface.  Once  one  accounts  for  the  syaaetry  of  the  free-space  Green's 
function  between  the  field  and  source  points.  It  becoaes  obvious  that  the 
Integration  over  the  surface  yields  an  expression  that  is  the  first  varia¬ 
tion  of  a  functional  J[p],  with  the  surface  velocity  held  constant.  The 
theorea  states  that  the  value  of  a  functional  J  that  depends  on  the  pressure 
field  is  stationary  to  Inflnlteslaal  variations  of  that  field  when  p  is  the 
true  pressure.  The  specific  foras  are 

f  J  -  0  (1) 


where 
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JtPl  -  IJJ  JJ  {k^[  n(O  n(|)  ]  p(f)  p(|) 

s  s 

-  [  n(f)  X  V^p(f)  ]  .  [  n(^)  X  v^p(0  G(M)  dA^  dA^ 

-  4*iwp  J  p(OU„(f)  dA^  (2) 

S 


Mote  that  the  variational  principle  govema  the  aurface  pressure  correspond 
Ing  to  a  specified  velocity  field,  so  the  noraal  velocity  v^  is  held 
constant  in  the  variation.  The  quantity  ^  functional  of  the  sur¬ 

face  velocity  whose  evaluation  entails  integrating  over  all  points  on  the 
surface,  specifically 


U„<C) 


pr|  JJ  v„(f>[ 


n(f) 


H-iA 

ar 


11 


gjO(f.l) 


dA, 


(3) 


In  these  expressions,  ar  is  the  distance  between  the  field  and  source  point, 
ar  •>  f  -I,  and  G(f,|)  denotes  the  nondimens ional  free-space  Green's  func¬ 
tion. 


G(M)  - 


-exp(ikar):  three -dimensions 
2^HQ^\kar):  two -dimens ions 


(A) 


Also,  PR  denotes  the  Cauchy  principle  value  of  the  integral,  which  is  ob¬ 
tained  by  excluding  from  the  region  of  integration  an  infinitesimal  patch 
centered  on  f 


Certain  features  of  the  expression  for  J[p],  Eq.  (2),  should  be  noted. 
All  gradients  of  pressure  at  any  point  on  the  surface  are  crossed  with  the 
nonsal  vector  at  that  point,  so  the  derivatives  represent  the  gradient  of 
the  pressure  tangentally  to  the  surface.  No  gradients  of  the  Green's  func¬ 
tion  appear  in  the  pressure  integrals,  but  the  normal  velocity  functional 
does  contain  such  a  gradient.  This  term  actually  has  the  same  sin¬ 
gularity  as  the  Green's  function,  because  the  dot  product  of  n(C)  and  f  -  | 
2  -  - 

behaves  as  r  as  ^  f.  As  a  conseq\ience ,  all  singularities  arising  in  the 
evaluation  of  G(f,0  are  integrable. 
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In  ganaral,  a  variational  principle  has  two  uses.  One  can  apply  the 
calculus  of  variations  to  obtain  the  governing  equations  for  the  systea. 

Such  an  eq\iation  is  known  as  the  Euler -Lagrange  equation  for  the  variational 
principle.  Not  surprisingly,  the  Euler -Lagrange  equation  obtained  from  Eqs. 
(1)  and  (2)  is  the  second  surface  integral  equation.  A  variational 
principle  may  also  be  used  as  the  basis  for  an  approximate  solution.  This 
is  the  way  in  which  it  was  used  in  the  project. 


II. 1  General  Implementation  of  SVP  for  Bodies  of  Bevolutlon 

The  latest  implementation  of  SVP  considers  an  arbitrary  excitation  of  a 
body  of  revolution.  The  shape  of  the  body  is  specified  by  describing  the 
cylindrical  coordinate  functions  of  its  shape  generator.  This  dependence  is 
expressed  in  parametric  form,  with  the  parameter  selected  as  the  arclength  s 
measured  along  the  generator  from  one  apex.  It  is  convenient  to  work  with 
nondimens ional  variables.  For  this  reason,  the  parameter  is  selected  as  a 
s/a,  where  the  reference  length  a  is  defined  as  the  largest  radius.  Thus, 
the  two  variables  locating  a  point  on  the  surface  are  the  nondimens lonal 
meridional  arclength  a  and  the  azimuthal  angle  f,  with  the  axial  and 
transverse  distance  for  this  point  given  by  specifying  aR(a)  and  az(a). 
Evaluation  of  the  SVP  integrals  requires  definition  of  the  second  surface 
point  f  which  is  specified  by  its  arclength  parameter  and  an  azimuthal 
angle  9  relative  to  the  first  point. 

The  arclength  and  azimuthal  angle  locating  a  point  on  the  surface  are 
also  used  to  describe  the  surface  pressure  and  surface  velocity.  The  2w  pe¬ 
riodicity  of  the  azimuthal  dependence  enables  one  to  expand  p  and  v^  in 
Fourier  series.  Using  the  complex  form  of  such  series  simplifies  some  of 
the  ensuing  steps,  so  the  pressure  and  normal  velocity  distributions  on  the 
surface  are  described  by 

p(f ,r)  -  I  pc^P(a,f)  exp(iwt)  +  C.C. 

v^(f ,r)  -  I  cT(o,^)  exp(i«t)  +  C.C.  (5) 
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where  C.C.  denotes  the  conplex  conjugate  of  preceding  terms,  and  the  spatial 
distributions  are 


-  5  ^  P^(o)  exp(ia«)  ^  I  X  ^m^“^  exp(im^) 


These  expressions  and  similar  ones  for  the  second  surface  point  are  sub¬ 
stituted  into  the  definition  of  J[p],  Eq.  (2).  This  enables  one  to  carry 
out  the  azisnithal  portion  of  each  surface  integral.  As  a  result  of  the 
axial  symmetry  of  the  shape,  the  azimuthal  harmonics  uncouple.  The  varia¬ 
tional  functional  is  thereby  found  to  have  the  form 


T  2  4 
J  p  c  a 


where  the  nondimens ional  functional  for  each  azimuthal  harmonic  is  given  by 

-  J  Iq  |ca<“./9)[kVR(o)R(^)R»(o)R'(/J)  -  m^z(o)z(^)]Pj^(o)P_^(^) 

+  D^(a,^)[[kVR(a)R(^)z'(a)z'(^)  -  m^R' (o)R' (/9)] P^(o)P_^(^) 

-  R(a)R(/»)P'(a)P'  (/9)1  -  E„(a,^)m[R' (o)R(^)P„(o)P'  (^) 

m  -m  J  m  *-  m  -m 

+  R(o)R' (^)P^(a)P  ^(i9)]|  da  d^  -  xW  R(o)T^(a)P  ^(a)  da 
2  £  t 

-  kV  F^(a,^)T^(a)P  da  dfi  (8) 


In  the  foregoing  a  prime  indicates  differentiation  with  respect  to  the  argu¬ 
ment,  and  i  -  max(s)/a. 


The  functions  C  (o,^),  D_(a,^),  E  (o,^),  and  F(a,/3)  are  integrals  over 
the  azimuthal  angle  of  the  Green's  function,  weighted  by  various  factors. 
Only  the  first  and  last  of  these  functions  need  to  be  computed.  The  others 
are  obtainable  from  recursive  identities. 


c  (o,^)  -  2  f  7  exp(ikar)  cos(m8)  dS 
■  Jo 


(9) 


D_(a.^) 


-I 


|[C^  ^(a.^)  +  C^i(o.i9)l  if  ■  >  0 


(10) 


C^(o.^)  If  ■  -  0 


E.(a,^) 


-1 


i[C^  ^(a.^)  -  C^i(a.^)]  If  m  > 
0  if  n  »  0 


(11) 


F^(o.^)  -  2R(a)R(^)  '(a)[R(a)  •  R(^)  cos  e] 

+  R'(a)[z(a)  -  z(^)]^  ^  exp(ikar)  cos(m8)  de  (12) 


where  the  distance  r  between  points  on  the  surface  Is  related  to  the  shape 
functions  by 


r  -  'Jr(o)^+  K{0)^  -  2R(o)R(^)cos  6  +  [z(a)  -  z(^)]^p/^  (13) 


Each  integrand  appearing  in  Eqs.  (9)  and  (12)  contains  a  1/r  singularity, 
whose  contribution  Is  isolated  as  elliptic  Integrals.  The  remainders  are 
coBiputed  by  conventional  numerical  methods. 


The  next  step  In  the  formulation  of  SVP  Is  to  express  the  vinknown  sur¬ 
face  pressure  as  Rltz  series  expansion,  whose  general  form  Is 


N 


B 


-I 


^mj 


(14) 


j-1 


idiere  the  coefficients  F  .  are  the  quantities  to  be  determined.  The  condl- 

mj 


tlons  imposed  on  the  N  basis  functions  are  minimal.  Necessary 


conditions  are  that  reconstruction  of  the  pressure  distribution  be  con¬ 
tinuous  everywhere,  including  the  boundaries  a  -  0  and  a  •  i,  which  are  the 
apexes  of  the  body  of  revolution.  Other  conditions,  such  as  symmetry  about 
any  of  the  coordinate  planes,  are  equivalent  to  natural  boundary  conditions 
in  structural  dynamics.  If  such  conditions  are  not  Inherent  to  the  basis 
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functions  selected  for  the  analysis,  they  will  nevertheless  be  exhibited  by 
the  SVF  solution. 


In  the  Initial  phases  of  the  project,  where  SVP  was  verified  against 
classical  analytical  results,  a  priori  knowledge  of  the  solution  was  usxially 
used  to  select  the  appropriate  basis  functions.  With  the  ultiaate  nerger  of 
SVP  with  structural  dynaaics,  such  knowledge  was  not  available.  The  wave- 
nuaber  based  foraulation  was  developed  to  address  this  sittiation.  Because 
the  arclength  aax(s)  is  the  extent  over  which  waves  aay  travel,  it  repre¬ 
sents  their  effective  aperture.  Correspondingly,  all  waves  Bust  be 
representable  as  half-range  Fourier  series  using  either  cosines  or  sines. 
Which  type  of  function  one  should  \ise  depends  on  the  aziauthal  haraonic  nua- 
ber.  For  a  -  0,  the  axisyaaetry  condition  at  the  apexes  is  satisfied  by 
cosines,  while  continuity  of  pressure  at  those  locations  requires  the  xisage 
of  sines  for  n  >  0.  Thus,  the  basis  functions  are  selected  as 


cos(jira)  if  a  -  0 
sin(jiro)  if  a  >  0 


(15) 


The  reason  for  referring  to  this  formulation  as  the  wave -number  based  SVP 
comes  from  replacing  the  trigonometric  functions  above  with  complex  exponen¬ 
tials,  and  then  combining  the  two  Fourier  series  for  meridional  and 
azimuthal  dependence  of  pressure.  The  form  of  this  combination  is 


CO  N 

■  I  I  e:cp[l[,rt  ±  -  „t]]  +  C.C.  (16) 

a— CO  j-1 


which  shows  that  the  P^j  coefficients  are  the  amplitudes  of  the  various 
presstire  waves  whose  rays  trace  spirals  along  the  wetted  surface.  Thus,  the 
results  obtained  from  this  formulation  of  SVP  are  the  analytical  analogs  of 
the  experimental  data  ultlsuitely  provided  by  acoustical  holography. 


The  final  step  in  the  formulation  of  SVP  involves  evaluating  the 
varioxis  functionals  J^,  Eq.  (8),  resulting  from  substitution  of  the  series 
expansion  for  surface  pressure,  Eq.  (14).  The  general  form  of  this  result 
is 
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N  N 

-III 

j-1  n-1 


2  2 

A  ,  P  ,  P  - 
mjn  Bj  -an 


N 

I 

j-1 


®  4  P  < 

■j  -“j 


(17) 


The  coefficients  A  .  and  B  .  are  obtained  from  double  integrals  expressing 
Bjn  Bj 

all  possible  position  pairs  along  a  Beridian.  Their  definition  is  the  saae 
as  the  corresponding  integrals  in  Eq.  (8)  above,  with  P^(a)  replaced  by 
etc.  Because  the  doBsin  of  Integration  for  both  sets  of  coeffi¬ 
cients  extends  over  the  singularity  at  a  the  doBain  is  broken  down  into 

squares  and  triangles  whose  sides  are  parallel  to  that  diagonal.  The 
nuBerical  values  are  obtained  by  application  of  a  Gaussian  rule  based  on 
nine  interior  points  for  each  subdoBain. 


Because  Eq.  (17)  describes  J  as  a  function  of  the  P  .  coefficients, 

a  BJ 

the  first  variation  is 


N 


I 


3J 


av 


J-1 


-mj 


5P  ,  -  0 
-mj 


(18) 


This  relation  must  hold  for  arbitrary  jP  ,  which  can  only  be  true  if  all 
partial  derivatives  vanish  Independently.  The  Batrix  form  of  these  deriva¬ 
tives  is 


(19) 


where  (P  }  contains  the  N  values  of  P  . . 

B  Bj 


Evaluation  of  the  coefficients  B  .  requires  a  description  of  the  Banner 

BJ 

in  which  the  surface  nonnal  velocity  depends  on  the  Beridional  arclength  for 
each  azimuthal  harmonic.  In  the  case  of  radiation  due  to  rigid  body  motion, 
such  dependence  is  readily  obtained  by  taking  the  component  of  the  rigid 
body-  velocity  in  the  direction  normal  to  each  point  on  the  surface. 
Evalxiation  of  radiation  from  elastic  bodies  requires  structural  dynamics 
considerations.  (Acoustic  scattering  is  treated  by  modifying  the  formula¬ 
tion  for  radiation  in  a  manner  that  will  be  discussed  later.) 
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11.2.  Structural  Dynauics  for  Platas  and  Shalls 


Many  aethods  for  analyzing  structural  notion  are  available,  including 
finite  eleaents,  finite  differences,  and  boundary  eleaents.  The  one 
selected  to  be  coid>ined  with  SVP  was  the  iwthod  of  assuaed  nodes,  which  aore 
properly  should  be  referred  to  as  a  Ritz  expansion.  This  aethod  is  closest 
to  SVP  froB  a  conceptual  standpoint,  because  both  represent  the  response 
over  the  entire  domain  in  terms  of  series  of  global  basis  functions.  (It 
would  not  be  difficult  to  modify  the  fomulation  to  accomodate  the  local 
basis  functions  associated  with  a  finite  element  description  of  the  struc¬ 
ture.)  The  deformation  of  a  plate  or  shell  is  described  by  displacement 
component  w  normal  to  its  aidsurface,  and  in-plane  displacements,  u  and  v. 
For  a  body  of  revolution  these  components  are  in  the  meridional  and 
azimuthal  directions.  Each  displacement  component  is  expanded  in  a  series 
comparable  to  that  used  to  describe  surface  pressure. 


•  N 

“  "  5  Z  Z  exp(iae) 

m— «  J-1 

•  N 

^  t  Z  Z  exp(ine) 

in—*  j-1 

«B  N 

w  -  i  ^  Z^  '^mj^*^^  ^mj^®^  exp(ime)  (20) 

m— «o  j-1 


where  the  various  are  admissible  basis  functions  for  the  displacement 
coiiq>onent  indicated  by  the  superscript.  (Admissible  functions  are  those 
that  satisfy  the  geometric  boundai^  conditions  and  have  the  appropriate  de¬ 
gree  of  differentiability.) 


The  equations  of  motion  governing  the  amplitudes  of  the  basis  functions 
are  obtained  by  using  the  series  to  construct  the  kinetic  and  strain  energy. 
The  virtual  work  must  also  be  described  in  terms  of  these  displacements, 
with  one  term  consisting  of  the  integral  of  (-p  Sv)  to  describe  the  effect 
of  the  acoustic  pressure.  This  procedure  leads  to  the  the  generalized 
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forces  associated  with  each  dlsplaceaent  conponent.  Hamilton's  extended 
principle  for  this  system  reduces  to  Lagrange's  equations  for  the  basis 
fxmctlon  amplitudes.  Redefining  these  amplitudes  by  factoring  out  the  har* 
BK>nlc  temporal  variation  leads  to 


(21) 


where  each  element  In  the  above  matrices  represents  an  N  dimensional  parti¬ 
tion.  The  various  K  and  M  represent  stiffness  and  mass  coefficients 
associated  with  the  displacement  components  Indicated  by  their  subscript, 
the  Q  terms  are  the  generalized  forces  due  to  the  dynamic  excitation,  and  F 
are  the  coupling  factors  describing  the  contribution  of  each  pressure  basis 
function  to  the  generalized  force. 


Because  the  normal  components  of  velocity  In  the  fluid  and  the  struc¬ 
ture  iBust  be  continuous  across  the  wetted  surface,  the  characterization  of 
the  normal  displacement  as  a  series  defines  the  terms  In  Eq  (6) ,  idilch 
Is  the  merlclonal  dependence  of  normal  velocity  In  the  a  th  azlonithal  har¬ 
monic.  Matching  the  meridional  dependencies  of  the  structure  and  fluid 
velocities  leads  to 

N 

j-1 

When  this  representation  of  the  meridional  variation  of  the  m  th  harmonic  of 
stirface  velocity  Is  used  to  evaluate  the  coefficients  In  Eq.  (17),  one 
finds 

N 

B  ,  -  y  A  .  W  (23) 

“j  L, 

j-1 
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wiiere  represent  coupling  factors  describing  the  napping  of  the  struc¬ 

tural  velocity  into  the  functional  space  twed  to  describe  the  pressure.  The 
result  of  using  the  foregoing  to  form  the  SVP  equatloiis  is 


(24) 


The  coupled  response  of  the  surface  pressure  and  the  notion  of  the  structure 
in  each  azinuthal  hamonic  nay  be  detemined  solving  the  SVP  equation  (24) 
sinultaneotisly  with  the  structural  dynanics  equation  (21) . 


II. 3.  Kodlficatlons  for  Scattering 

The  developnent  thus  far  addressed  radiation  problens.  A  simple 
nodification  extends  the  treatment  to  cases  %fhere  the  excitation  is  an  inci¬ 
dent  acoustic  signal.  The  key  is  to  decompose  the  pressure  into  incident 
and  scattered  parts.  Then  the  scattered  pressure  on  the  surface  satisfies 
the  sane  set  of  surface  integral  eq\iations  as  the  total  surface  pressure  in 
a  radiation  problem,  provided  that  the  normal  velocity  is  defined  as  the 
particle  velocity  in  the  scattered  signal.  Continuity  of  normal  velocity  on 
the  surface  requires  that  this  particle  velocity  is  the  difference  between 
the  structure's  velocity  and  the  nomal  component  of  particle  velocity  in 
the  Incident  wave.  Euler's  equation  gives  the  latter,  so  that 

n-Vp  -kaw-n’Vp  (25) 

It  follows  that  replacing  any  term  traceable  to  n  ■  Vp  in  the  SVP  equations 
by  the  rigjht  side  of  the  above  leads  to  relations  governing  the  scattered 
part  of  the  pressure  field. 

In  order  to  apply  this  result,  the  incident  pressure  and  particle 
velocity  are  expanded  in  azimuthal  Fourier  series. 

« 

inc.r.  1  2  Y""  „inc,  .  r., 

P  (D  -  2  exp(ime) 


(26) 


n(f)  •  Vj,p  (f)  -  j  />c 


I 


T^"®(a)  •xpdne) 


This  converts  the  SVP  eqiistlons  to 


[  *.  ]lf.l  -  A.  -  kV(l>*“) 


2  2.„lnc, 


inc 

idiere  (B  )  represents  the  values  of  B  .  corresponding  to  a  surface 
^  inc 

velocity  irtiose  meridional  dependence  is  T  (a).  In  addition,  the  presence 

B 

of  an  incident  pressure  field  adds  to  the  generalized  forces  driving  the 
normal  displacement.  This  leads  to  an  additional  generalized  force,  which 
is  described  by  replacing  Qy  solution  of  the  SVP  and 

structural  dynamics  equations  in  this  approach  yield  the  displacement 
amplitudes  and  the  scattered  part  of  the  pressure  field.  The  total  surface 
pressure  may  then  be  determined  by  superposing  the  incident  field's  con¬ 
tribution  to  each  azintuthal  harmonic. 


III.  RESBARCH  TOPICS  AND  RESULTS 

The  following  discusses  the  themes  of  the  various  research  phases  of 
the  project.  The  numbers  appearing  at  the  end  of  each  section  refer  to  the 
papers  listed  in  the  next  section. 

III. I.  Radiation  from  Vibrating  Rigid  Bodies  (papers  #  1-8) 

Each  research  phase  in  which  the  SVP  approach  was  generalized  began 
with  the  task  of  evaluating  the  surface  pressure  on  a  body  whose  motion  is 
prescribed.  The  first  cases  to  be  pursued  were  a  transversely  oscillating 
thin  rigid  disk,  a  sphere  whose  surface  is  executing  a  breathing  mode  vibra 
tion,  and  a  rigid  sphere  and  a  flat-ended  cylinder  oscillating  parallel  to 
the  axis  of  symmetry.  The  basis  functions  for  pressure  were  initially  sug¬ 
gested  by  the  available  analytical  solutions,  rather  than  the  sinusoidal 


-  Page  15  - 


fonu  associated  with  the  wavenusber  foraulation  of  SVP.  A  study  of  the  os¬ 
cillating  rigid  disk  was  also  carried  out  using  eleaental  basis  functions 
coaiparable  to  those  enployed  for  boundary  elenents. 

Each  of  the  foregoing  cases  constitutes  an  axisyaoietric  systea.  The 
extension  of  SVP  to  nonsymetric  cases  addressed  a  heai-capped  cylinder,  and 
a  spheroid  undergoing  rigid  body  aotion  in  the  axial  and  beaawise  direction, 
as  well  as  a  rotational  oscillation  about  the  centroid.  The  foraer  problea 
provided  an  iaportant  consistency  check.  If  one  sets  the  length  of  the 
cylindrical  section  to  zero,  the  systea  is  a  sphere.  The  pressure  distribu¬ 
tion  for  beaawise  translation  is  not  axisyaaetric ,  whereas  axial  translation 
represents  an  axisyaaetric  noraal  velocity  distribution.  It  was  verified 
that  the  SVP  results  for  the  beaawise  case  gave  the  saae  surface  pressure  as 
that  obtained  for  axial  aotion,  and  both  results  were  verified  to  give  the 
saae  result  as  the  analytical  solution  for  the  axial  translation.  It  was  in 
this  phase  of  the  research  program  that  the  wavenun^er  foraulation  of  SVP 
was  developed. 

Addressing  problems  of  this  type  was  very  useful  for  a  number  of 
reasons.  Because  structural  dynamics  equations  were  not  an  issue,  it  was 
easier  to  debug  the  implementation  of  the  SVP  equations.  Furthermore,  rigid 
body  notions  excite  only  the  m  -  0  and  a  -  1  azimuthal  harmonics,  so  it  was 
not  necessary  to  worry  about  convergence  and  nxiaerical  error  associated  with 
high  haraonics.  Another  aid  to  validating  the  analysis  was  the  fact  that 
analytical  solutions  were  available  in  some  cases.  In  addition,  studies  of 
slender  cylinders  and  prolate  spheroids  led  to  some  interesting  results 
regarding  the  the  effects  of  curvature  on  a  surface.  In  every  case,  conver¬ 
gence  of  the  pressure  amplitudes  led  to  results  that  differed  by  less  than 
IX  from  those  of  analytical  solutions.  Similar  close  agreement  for  a  heai- 
capped  cylinder  was  obtained  from  a  comparison  with  the  predictions  of  the 
CHIEF  boundary  element  program.  It  is  important  to  realize  that  convergence 
of  the  series  occurs  only  if  the  higher  order  terms  in  series  are  very 
small.  It  follows  that  one  can  examine  the  magnitude  of  the  hi^er  order 
terms  to  determine  that  the  solution  has  converged- - there  is  no  need  to  per¬ 
form  an  analysis  using  a  diiferent  series  length  to  determine  the  adequacy 
of  a  result. 
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III. 2.  Eadlatloa  froa  Blastie  BodlM  (papers  #  9-22) 

The  research  plan  was  to  use  the  rigid  body  notion  problea  to  general¬ 
ise  the  SVP  foraulation.  Once  each  effort  was  finished,  the  next  task  was 
to  iapleaent  the  structural  dynanics  aspects  for  that  version.  Thus,  the 
first  case  where  structural  dynanics  was  considered  was  a  circular  elastic 
nenbrane  or  flexural  disk  whose  boundary  is  supported  by  an  annular  baffle. 
The  width  of  this  baffle  could  vary  fron  xero  to  infinity.  Only  axisynnetric 
excitations  were  considered.  The  basis  functions  for  the  structure  in  these 
studies  were  selected  as  the  In- vacuo  nomal  nodes.  Several  types  of  basis 
functions  for  the  surface  pressure  representation  were  addressed.  It  was 
shown  that  the  basis  functions  used  to  study  the  rigid  body  case  lead  to 
ill-conditioning  when  a  large  number  of  functions  are  retained,  because  the 
higli  order  functions  are  not  linearly  independent  from  a  computational  view¬ 
point.  Bessel  and  sinusoidal  basis  functions  were  shown  to  agree  well  with 
the  analytical  solution  derived  by  Alper  and  Magrab,  as  were  boundary 
element-type  functions  that  spanned  small  sxibdomains.  A  further  study  of 
the  elastic  plate  considered  the  effect  of  annular  stiffeners.  This  problem 
introduced  the  concept  of  component  mode  analysis,  which  was  subsequently 
implemented  by  Igusa  and  Acheiibach  of  Northwestern  University  as  a  fundamen¬ 
tal  aspect  of  their  ONR  sponsored  research  program. 

The  next  elastic  structure  to  be  addressed  was  a  spheroidal  shell  of 
arbitrary  aspect  ratio.  Two  formulations  of  the  fluid -structure  coupling 
were  derived.  The  first,  which  was  termed  the  "direct  method"  followed  the 
formulation  described  in  the  previous  section.  The  second  was  called  the 
"nodal  method,"  because  the  assumed  modes  method  is  used  to  evaluate  the  in- 
VACUO  natural  frequencies  and  mode  shapes.  These  free  vibration  properties 
were  then  used  the  basis  functions  for  the  SVP  formulation.  The  results 
for  free  vibration  were  shown  to  be  extremely  accurate,  and  to  be  more  com¬ 
prehensive  than  any  previ^'^isly  available.  An  interesting  feature  was  shorn 
to  be  a  coelescence  of  eigenvalues  at  certain  aspect  ratios,  which  launched 
a  separate  research  phase  that  will  be  discussed  later.  The  only  data 
available  for  validation  in  the  fluid- loaded  problem  was  for  a  spherical 
shell.  The  SVP  results  were  shown  to  be  very  close  to  Hayek's  analytical 
solution. 
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An  •lutic  plat*  in  an  annular  baffla  was  the  first  study  to  renove  the 
restriction  to  axisyaaetric  situations.  The  plate  was  s\ibjected  to  an  ec¬ 
centric  concentrated  hanonic  force  nonal  to  the  surface.  This  study  was 
the  first  to  inplMMnt  the  Fourier  series  representation  of  the  aziauthal 
dependence.  No  prior  work  had  considered  this  problen,  but  the  results  were 
examined  for  self-consistency  and  physical  reasonableness.  At  the  present 
time,  the  corresponding  SVP  iq>lenenCation  for  shells  is  under  development. 

111. 3.  Scattering  from  Rigid  Bodies  (paper  #  23) 

It  was  decided  to  delay  considering  acoustic  scattering  \mtil  the  SVP 
formulation  was  generalized  to  handle  non-axisymaetric  problems.  The  inci¬ 
dent  signal  in  each  SVP  analysis  was  a  plane  acoustic  wave  at  arbitrary 
incidence.  Results  have  been  obtained  for  scattering  from  flat  and  heml- 
capped  cylinders,  and  from  prolate  spheroids  whose  surface  is  stationary. 
These  were  compared  to  numerical  solutions  obtained  from  the  CHIEF  and  non- 
symmetric  SHIP  programs.  In  every  case  the  results  were  virtually  identical 
to  the  nuBwrical  solutions.  A  consistency  check  was  also  carried  out  for 
the  case  of  a  sphere,  where  the  solution  should  yield  the  same  pressure  dis¬ 
tribution  as  that  predicted  analytically,  independent  of  the  angle  of 
incidence . 

111. 4.  Ihm  Interior  Cavity  Problem  (papers  #  21,  22,  24,  25) 

Because  the  stationary  principle  for  an  SVP  formulation  is  derived  from 
the  second  surface  integral  equation,  there  is  no  unique  solution  at  fre¬ 
quencies  that  match  the  eigenfrequencies  for  the  interior  Neumann  problem 
(zero  velocity  on  the  boundary) .  A  study  was  performed  to  assess  whether 
the  over-determined  formulation  employed  in  CHIEF  could  be  used.  This  en¬ 
tailed  adding  auxiliary  constraint  equations  requiring  that  the  interior 
pressure  vanish  at  selected  points  or  on  average  over  an  interior  surface. 
Although  the  concept  worked  well,  it  has  not  been  implemented  into  the 
general  SVP  formulation.  This  decision  stemmed  from  observations  in  the 
analysis  of  radiation  from  a  spheroidal  shell.  That  work  showed  that  even 
if  one  were  to  come  very  close  to  one  of  these  frequencies,  the  effect  would 
be  either  very  minor  or  imperceptible.  The  apparent  reason  for  this  result 
is  the  limited  size  of  an  SVP  model,  i.e.  the  relatively  small  nuid>er  of 
terms  required  to  obtain  an  accurate  solution.  As  a  consequence,  the  system 
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of  oiiuotions  roaolns  well -conditioned  even  if  one  ia  within  O.OIX  of  the 
spurious  resonence  frequency.  This  is  an  iaportsnt  factor  because  the 
decoaposition  of  the  problem  into  azimuthal  harswnics  means  that  the 
spurious  resonance  frequencies  will  occur  in  relatively  equal  density  at 
each  azimuthal  harmonic.  In  contrast,  conventional  botmdary  element  for¬ 
mulations  encounter  an  increasing  density  of  spurious  frequencies  as  the 
excitation  frequency  increases,  because  they  treat  all  azimuthal  harmonics 
simultaneously. 

III. 5.  ligenvmlus  Veering  Ihenomens  (papers  #  26-33) 

The  analysis  of  the  in-vacuo  modes  of  a  spheroidal  shell  disclosed  that 
at  certain  length  to  radius  ratios,  the  natural  frequencies  of  the  flexxural 
and  extensional  branches  become  very  close.  Further  Increase  of  this  aspect 
ratio  led  to  modal  vibration  patterns  that  exhibited  localized  displacement 
patterns.  Comparable  phenomena,  which  are  grouped  under  the  term  of 
"eigenvalvw  veering  phenomena",  had  been  studied  for  other  systems. 

However,  all  such  work  failed  to  describe  the  actual  extreme  sensitivity  of 
the  mode  shapes  as  the  system  parameter  is  varied  through  the  zone  where 
eigenvalue  veering  occurs.  A  general  treatment  was  derived  as  part  of  the 
research  program,  and  then  applied  to  the  spheroidal  shell  problem.  It  was 
shown  that  the  localized  modes  cannot  be  omitted  from  the  modal  version  of 
SVP,  even  if  their  natural  frequencies  are  much  higher  than  the  excitation 
frequency.  This  result  apparently  stems  from  the  extreme  efficiency  with 
which  these  modes  radiate. 

Another  study  of  the  effect  of  eigenvalue  veering  phenomena  on  acoustic 
radiation  was  carried  out  using  a  two-dimensional  flat  plate  in  an  infinite 
baffle.  Addressing  this  situation  necessitated  extending  SVP  to  handle  two- 
dimensional  (plane  strain)  problems.  That  formulation  was  validated  by 
comparing  the  SVP  results  for  a  single -span  plate  to  the  integral  transform 
solution.  Then  the  formulation  was  applied  to  a  two-span  plate  %fhose  in¬ 
dividual  spans  are  nearly  uncoupled  by  a  stiff  torsional  spring  acting  at 
the  center  support.  The  in-vacuo  analog  of  this  problem  had  served  as  a 
protoQrpe  for  more  than  a  decade  for  studies  of  eigenvalue  veering.  The 
startling  result  that  emerged  form  this  work  was  that  a  shifting  the  posi¬ 
tion  of  the  center  support  by  as  little  as  IX  of  the  span  length  can  raise 
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th«  total  rodlatod  power  by  a  factor  of  two.  Because  a  coapllcated  systea 
stich  as  aubaarine  Is  sure  to  have  natural  frequencies  that  are  very  close, 
It  Is  reasonable  to  anticipate  that  conparable  eigenvalue  veering  phenoaena 
will  occtir  in  those  systeas.  This  has  significant  iaplications  for  con¬ 
structing  aatheaatical  aodels,  because  it  suggests  that  little  details  of 
construction  aay  have  a  strong  influence  on  acoustical  perforaance. 


IV.  CLOSnE 

This  research  project  has  resulted  in  developaent  of  the  SVP  technique 
for  aodeling  the  coupled  acoustic-structure  response  of  stibaerged  bodies. 

The  technique  requires  detemination  of  a  eoaparatively  saall  nuBA>er  of 
variables.  In  the  wavenuaber  version  of  SVP,  these  variable  are  par¬ 
ticularly  aeaningful  froa  a  physical  standpoint,  becaxise  they  represent  the 
aaplltudes  of  helical-type  waves  that  propagate  through  the  structure  and 
through  the  fluid.  SVP  solutions  have  the  i8q>ortant  property  of  allowing 
one  to  ascertain  that  the  results  are  convergent  and  represent  the  response 
to  a  satisfactory  level,  without  requiring  evaluation  of  a  smre  refined 
aodel  for  coaparison.  A  nuaber  of  generalizations  ultiaately  yielded  a 
technique  that  is  equally  suitable  for  radiation  and  scattering.  The  only 
capability  requiring  further  developaent  is  analysis  of  non- axisyaae trie 
probleas  in  curved  shell  structures.  Because  of  the  relatively  saall  nuaber 
of  unknowns  to  be  deterained,  and  the  fact  that  those  unknowns  are  physi¬ 
cally  aeaningful,  the  technique  provides  physical  insist  into  the  response 
of  siibaerged  structures.  These  features  aake  SVP  very  useful  as  a  research 
tool,  as  well  as  a  tool  for  design  paraaeter  studies. 
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